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Charge and spin density distributions are studied within a nano-ring structure endowed with 
Rashba and Dresselhaus spin orbit coupling (SOI). For a small number of interacting electrons, 
in the presence of an external magnetic field, the energy spectrum of the system is calculated 
through an exact numerical diagonalization procedure. The eigenstates thus determined are used to 
estimate the charge and spin densities around the ring. We find that when more than two electrons 
are considered, the charge-density deformations induced by SOI are dramatically flattened by the 
Coulomb repulsion, while the spin density ones are amplified. 
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I. INTRODUCTION 

The possibility of controlling the flow of the electron 
spins in semiconductor structures by external electric 
means through spin-orbit interaction (SOI) has domi- 
nated the recent past of spintronics research. This fun- 
damental principle, first explored in the Datta-Das spin 
transistor configuration^ has been guiding a sustained 
effort in understanding all the phenomcnological impli- 
cations of this interactions on systems of electrons. The 
coupling between spin and orbital motion results either 
from the two-dimensional confinement (Rashba)^ or from 
the inversion asymmetry of the bulk crystal structure 
(Dresselhaus)^ The usual expression of the the spin-orbit 
Hamiltonian Hso retains only the linear terms in the 
electron momentum p = (p x ,Py) and is given by 

a j3 

H S = -ziPxPy ~ VyPx) + jKPxPx ~ °yPy) ■ (1) 

The Rashba and Dresselhaus coupling constants are a 
and /3, respectively, while <J x ,y,z are the Pauli matrices. In 
general, the two interactions are simultaneously present 
and often have comparable strengths. While a, the cou- 
pling constant of the Rashba interaction, can be modi- 
fied by external electric fields induced by external gates, 
the strength of the Dresselhaus SOI, (3, is fixed by the 
crystal structure and by the thickness of the quasi two- 
dimensional electron system^^ In many situations of in- 
terest, an additional energy scale is introduced by the 
Zeeman interaction of the electron spin with an external 
magnetic field, proportional to the effective gyromagnetic 
factor, g* , which depends on the material energy-band 
structure. While g* = —0.44 is very small in GaAs, it 
can be more that 100 times larger in InSb. 

The interplay between the two types of SOI, which 
have competing effects on the precession of the electron 
spin as they rotate it in opposite directions, and the 
Zeeman splitting, which minimizes the energy by align- 
ing the spin parallel to the external field, determines 
the ground state polarization of the electron system and 



the characteristics of spin transport. The investigation 
of such problems in mesoscopic rings has been pursued 
intensively by several authors^— In particular, it was 
noticed that, in the absence of the Coulomb interac- 
tion among the electrons, the interference between the 
Rashba and Dresselhaus precessions relative to the or- 
bital motion, leads to the creation of an inhomogeneous 
spin and charge distribution around the ringj£ The charge 
inhomogeneity created in this situation has a symmetric 
structure with two maxima and two minima around the 
ring and will be called here a charge-density deformation 
(CDD). The effect of the Coulomb interaction on this 
type of charge distribution has been considered for two 
electrons. It was obtained that, on account of the elec- 
trostatic repulsion, the two electrons become even more 
localized in the potential minima associated to the CDD, 
leading to an amplitude increase 

In this work we obtain an estimate of the effect of the 
Coulomb interaction on the charge and spin distribution 
associated with N = 2,3 and 4 electrons in a ring with 
SOI coupling by an exact diagonalization procedure that 
uses the configuration interaction method. Our results 
indicate that when the number of electrons increases, the 
mutual repulsion leads to more uniform charge distribu- 
tion around the ring, generating a dramatically flattened 
CDD. In contrast, the spin-density deformation (SDD) is 
amplified by the Coulomb effects. This can be explained 
by the appearance of a stronger repulsion between same 
spin electrons, leading to more favorable spin orienta- 
tions. 



II. THE RING MODEL 

The system of interest in our problem is a two- 
dimensional quantum ring of exterior and interior radii, 
R ex t and Ri n t respectively. The ring is placed in a per- 
pendicular magnetic field B associated in the symmetric 
gauge with a vector potential A = B/2(—y,x,0). The 
single-particle Hamiltonian of an electron of momentum 
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FIG. 1: (Color online) The discretized ring with Ri„t — 
0.8R ex t, and 10 radial x 50 angular sites. The sites are shown 
with circular points. The thin dotted lines connection sites 
are for guiding the eye. 



p = —ihV + eA and effective mass m* is written as the 
sum of an orbital term Ho = p 2 /2m*, a Zeeman contri- 
bution Hz = (l/2)<?* hbBct z and the spin-orbit coupling 
given in Eq. (pj. The ensuing expression, 



H = Hq + H z + H, 



so 



(2) 



is discretized in a standard manner— ill on a gricU^ 
defined by N r radial and N v angular sites, as shown 
in Fig. [TJ The radial coordinate of each site is = 
Rext — (k — l)5r, with k = 1, 2, ...,N r , while Sr = 
(Rext — Rint)/(N r — 1) is the distance between adja- 
cent sites with the same angle. Similarly, the angular 
coordinate is (pj = (j — l)Sip, where j = 1,2, ...,N V and 
5ip = 2ir/N v is the angle between consecutive sites with 
the same radius. The Hilbert space is spanned by the 
ket-vectors \kja), where the first integer, k, stands for 
the radial coordinate, the second one, j, for the angular 
coordinate, and a = ±1 denotes the spin projection in 
the z direction. 

In this basis {\kja}}, the matrix elements of the orbital 
Hamiltonian arc given by: 



(kja\H \k ! j'a') = 



TS n 



tin ~f~~ H - f> ** 



!,2 ( r k 



\4R e 



dkk-Sjj' 



(3) 



tip + ts 



Skk'Sjj'+i + t r 5kk'+iS 




T = h 2 /(2m*Rl xt ) is the energy unit, while R ext is the 
length unit. In T units, we obtain t v = [R ex t/ '{fk^'p)] 2 
the angular hopping energy, t r = (R ex t/Sr) 2 the radial 
hopping energy, and tg = hcB / (m*T) the magnetic cy- 
clotron energy (h.c. denotes the Hcrmitian conjugate.) 



In the same basis, the Zeeman Hamiltonian is simply 
diagonal in the spatial coordinates, 

(kja\H z \k'fa') = ^Tt Bl {a z ) aa ,5kk>5 or , (4) 

where 7 = g*m* /(2m e ) is the ratio between the Zeeman 
gap and the cyclotron energy, m e being the free electron 
mass. 

For the spin-orbit Hamiltonian we obtain: 



(kja\Hso\k'j'cr') = ^Tt a 
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(5) 



where t a = a/(R ex tT) and tp = (3/(R ext T) are the 
two types of spin-orbit relative energies, while cr r (p) — 
a x cos tp + <iy sin ip and u v (p) = —<J X sin ip + a y cos (p are 
the radial and angular Pauli matrices, respectively. We 
used the slightly shorter notations a 3 r = a r (pj) and 
a 3 = a v (ipj) for the matrices at the particular angles 
on our lattice. The Rashba spin-orbit terms are all in- 
cluded in the first square bracket. The Dresselhaus terms 
are very similar to the Rashba ones, being given by the 
substitutions indicated in the second square bracket. 

The single particle states of the noninteracting Hamil- 
tonian ©j Hip a = taipa, are computed as eigenvalues 
and eigenvectors of the matrices (0)-©, tp a (rk,<Pj) = 
So- ^ a,a{k, j)\a) . where ^/ a ,a(k,j) are c— numbers. 

In the basis provided by {tpa} the interacting many- 
body Hamiltonian is written in the second quantization 
as 



^ Vabcdc\c\c d C c , 



(0) 



abed 



where c\\ and c a are the creation and annihilation opera- 
tors on the single-particle state a. The matrix elements 
of the Coulomb potential V(r — r') = e 2 /(K|r — r'Q, k be- 
ing the dielectric constant of the material, are in general 
give by 

V abcd = (V> (r)^b(r')|V(r - r')|Vc(r)Vd(r')) ■ (7) 

In the present discrete model the double scalar product 
is in fact a double summation over all the lattice sites 
and spin labels 



V abcd = Tt c £ KAkj)^i a ,(k'f ) 



R, 



kja 
k'j'a' 



r,fc - r fk , 



(8) 
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The new energy parameter introduced by the Coulomb 
repulsion is tc = G 2 /{KR ext T). In the above summation 
over the sites, the contact terms (k = k',j = j') are 
avoided, as their contribution vanishes in the continuous 
limit. 

The many-body states $^ are found by solving the 
eigenvalue problem for the Hamiltonian ([5]), 

H$ M = . 

A potential solution of the equation is written in the con- 
figuration interaction representation^— as a linear com- 
bination of the non-interacting system eigenstates (Slater 
determinants), 

$ M = ^ c ^|a)> (9) 

a 

with {|a} = \i°, i%, ig-)} where i" = 0, 1 is the occu- 
pation number of the single-particle state tp a and K is 
the number of single-particle states considered. The oc- 
cupation numbers are listed in the increasing energy 
order, so ck is the highest energy of the single-particle 
state included in the many-body basis. For any \a) we 
have i™ = N, which is the number of electrons in the 
ring. It is straightforward to derive the matrix elements 
of l-l aa ' using the action of the creation and annihila- 
tion operators on the many-body basis. In practice Eq. 
@ is convergent with K for a sufficiently small number 
of electrons, and sufficiently small ratio of Coulomb to 
confinement energy, tc- This procedure, also known as 
" exact diagonalization" , does not rely on any mean field 
description of the Coulomb effects, like Hartrcc, Hartrce- 
Fock, or DFTJ£ 

To be able to carry the numerical calculations in a rea- 
sonable amount of time, we choose a small ring of radii 
Rext = 50 nm and Ri n t = 0-8R ex t, containing N < 4 
electrons. The discretization grid has 10 radial and 50 an- 
gular points (500 sites) , as shown in Fig. [TJ Two common 
semiconductor materials used in the experimental spin- 
tronics are used for the selection of the material constants 
needed: In As with m* = 0.023?7i e , g* = -14.9, k = 14.6, 
and estimated (or possible) values for the spin-orbit in- 
teractions a ~ 20 and (3 ks 3 meVnm; InSb with m* = 
0.014m e ,3* = -51.6, k = 17.9, and a « 50, /3 30 
meVnm4^ The relative energies which we defined are: 
for InAs t a = 0.60, tp = 0.09, t c = 2.9,7 = -0-17; for 
InSb t a = 0.92, tp = 0.55, t c = 1.5,7 = -0-36. In 
our calculations we have considered material parameters 
somewhere in between these two sets: t a = 0.7, tp = 
0.3, t c = 2.2,7= -0.2 

III. RESULTS 
A. Single particle calculations 

In the absence of the SOI, (a = (3 = 0), the single 
particle Hamiltonian ^ shares its eigenstates ijj a with 
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FIG. 2: (Color online) (a) The lowest 12 energies of the single 
particle states vs. the magnetic energy ts- The solid (red) 
and the dashed (green) lines show the states with positive 
and negative parity, respectively, (b) The expected value of 
the spin projection along the z direction S z , in units of h/2, 
for the first four states on the energy scale, (c) The standard 
deviation A c of the charge distribution around the circle with 
radial site index k = 6, for the first four energy states: ground 
state (GS), first, second, and third exited states (ESI, ES2, 
ES3). The same association of line types with states is used 
in panel (b). 



the z components of the angular momentum L z and spin 
S z = h<r z /2. In the presence of only one type of SOI, ei- 
ther a ^ or j3 ^ 0, the Hamiltonian commutes with the 
z component of the total angular momentum, L z + S z , 
which is conserved. When both a ^ and (3 ^ 0, 
the angular momentum is no longer conserved. How- 
ever, tp a continue to be eigenstates of the parity operator 
V = Tla z , n being the (three dimensional) spatial in- 
version operator. Indeed, the general Hamiltonian ([2]) 
commutes with V, which can be easily verified by using 
np = — pn and the commutation rules of the Pauli ma- 
trices. So in general Vip a = sip a , and thus the parity 
,s = ±1 of any state a is conserved, i. e. it is indepen- 
dent on the magnetic field. In particular, when a = /3 
and g* = 0, all states become parity-degenerate at any 
magnetic field£ ' 10 ' 17 We identify the parity of the single 
particle states calculated on our discrete ring model by 
looking at the relation tp a ^ a (k,j) = saTp a ,a(k,j)) where 
(k,j) and (k,j) arc diametrically opposed sites, with an- 
gular coordinates ipj — (fj + ir. 

In Fig. [2] we show the single particle states energy for 
< ts < 10 (units of T), which corresponds to a mag- 
netic field strength between and 1.32 Tesla. Further 
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increment of the magnetic field requires an augmentation 
of the number of sites on the ring in order to maintain the 
discrete model as a reasonable approximation of a physi- 
cally continuous ring. At zero magnetic field all states are 
parity degenerate, which is just the ordinary spin degen- 
eracy. The degeneracy is in general lifted by a finite mag- 
netic field. There are, however, some particular values of 
where the degeneracy persist. This situation is repre- 
sented in Fig. [2ja) by all intersection points of two lines 
corresponding to the two possible parities. Such intersec- 
tions do not occur between states with the same parity. 
Due to the spin-orbit coupling, the orbital momentum 
of one state depends on the spin of the other state and 
vice versa, and, consequently, states of same parity do 
in fact interact and thus avoid intersections.— Although 
in Fig. [U^a) many states represented by the same line 
type apparently cross each other, in reality there are al- 
ways tiny gaps between them, similar to those visible at 
ts ~ 2 between the first and the second excited states 
or at ts ~ 5.5 between the first, second, and third ex- 
cited states. The magnitude of these gaps depends on the 
g— factor, reducing in size for a smaller g* parameter. 

In Fig. [5£b) the evolution of the expected spin in the 
z direction, S z = h(i^ a \o' z \ipa) /2 is presented for the first 
four states in the energy order. One can see how the 
spin flips for states avoiding the crossing, like those with 
negative parity at ts ~ 2 (dashed and dotted lines in 
Fig. [Ib)). 

Only one type of SOI, either Rashba or Dresselhaus, 
is sufficient to avoid the crossing of states with the same 
parity, but in this case the charge and the spin densi- 
ties are uniform around the ring. When both SOI types 
are present the charge and spin densities become nonuni- 
form. This situation is equivalent with the presence of 
a potential with two maxima and two minima around 
the ring, having reflection symmetry relative to the axes 
y = x (or ip = 7r/4, corresponding to the crystal direction 
[110]) and y = —x (or ip = — 7r/4, corresponding to the 
crystal direction [110]) The amplitude of the CDD 
is illustrated in Fig. [2jc) where the standard deviation 
(in the statistical sense) A c of the charge density cal- 
culated around one circle on the ring, close to the mean 
radius, with radial site index k = 6, is plotted for the low- 
est four energy states. The density deformation occurs 
on account of the two combined SOI types which lead 
to spin interference and additional interaction between 
states with the same parity. Consequently, the ampli- 
tude of the CDD for a certain state is maximum at those 
magnetic fields where the parity degeneracy is lifted (the 
state avoids a crossing with another state of the same 
parity). In Fig.[5Jc) this is clearly seen at ts ~ 2, for the 
excited state. In this example the CDD in the ground 
state is very weak. The sharp peak at ts ~ 0.4 indicates 
the existence of a narrow gap between the 4-th and 5-th 
energy states that avoid crossing. 




FIG. 3: (Color online) Energy spectra of the first 12 states 
for N = 2, 3, and 4 electrons without Coulomb interaction, 
to — in panels (a),(c),(e), and with Coulomb interaction, 
to = 2.2, in panels (b),(d),(f). The solid (red) and the dashed 
(green) lines show the states with positive and negative parity, 
respectively. 



B. Many particle calculations 

In the following considerations, we will include more 
than one electron. In Fig. [3]we compare the energy spec- 
tra for the first 12 states vs. the magnetic energy for 
N = 2,3, and 4 electrons, without and with Coulomb 
interaction. Since the Coulomb interaction is invariant 
at spatial inversion (and independent on spin) the par- 
ity s is also conserved in the many-body states. Spectra 
drawn for tc = and tc = 2.2 have similar features. 
The interacting spectrum presents a shift to higher ener- 
gies, on account of the additional Coulomb energy, and a 
slight increase of the gaps at high magnetic fields. More- 
over, the crossings and the anti-crossings (points where 
the crossings were avoided) of the energy levels have a 
tendency to shift slightly to higher magnetic fields. Sim- 
ilarly, the gap between the ground state and the excited 
states increases at high ts- 

The total spin S z for each of the first three energy 
states is shown in Fig. 01 At zero magnetic field, for an 
even number of electrons, here N = 2 or N = 4, the 
ground state is non-degenerate and has total spin S z = 
0, i. e. the spin-up and spin-down states of individual 
electrons compensate. When the field is applied, the first 
spin flip in the interacting ground state, as well as the 
spin saturation, occur at lower magnetic fields than in 
the absence of the Coulomb repulsion. This is a result of 
the mixing of spin states with the same parity produced 
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FIG. 4: (Color online) The total spin projection in the z 
direction, in units of h/2, for the many body states with JV = 
2, 3, 4 electrons. Without interaction, i. e. to — in panels 

(a) ,(c),(e), and with interaction, with tc = 2.2, in panels 

(b) ,(d),(f). Only the first three states are shown here, the 
ground state (GS), the 1-st excited states (ESI), and the 2-nd 
excited state (ES2). The magnetic energy ts varies between 
and 10 and the lines showing the excited states are shifted 
to the right, for clarity. 

by the interaction. For N = 3 the ground state is spin 
(double) degenerate at zero field. In the presence of the 
Coulomb interaction, the total spin in the ground state 
and the higher state is reversed. 

Similar to the case of one electron (N = 1, Fig. [5]), the 
charge deformation of each state is maximized for those 
magnetic fields where the state has an anti-crossing (or 
repulsion) with another state of the same parity. The 
charge deformation parameter A c is shown in Fig. [5] For 
N = 2 the amplitude of the CDD increases with the 
Coulomb interaction. There is a simple reason for that: 
the potential associated with the charge deformation has 
two minima diametrically opposite on the ring and each 
of the two electrons tends be localized in one of these min- 
ima. The mutual Coulomb repulsion fixes the electrons in 
those places better The situation changes, however, 
for TV > 2. The Coulomb forces spread the electrons 
differently, more or less uniformly, such that the charge 
deformation created by the SOI is drastically reduced. In 
other words, the associated potential is strongly screened 
even by one extra electron above TV = 2. This effect can 
be clearly seen in Fig. [5l comparing panels (c) with (d) 
and (e) with (f ) . The vertical scale of panels (d) and (e) 
has been magnified three times, for visibility. 

In principle, the screening of the charge deformation 
is not particularly related to the spin-orbit effects. SOI 



FIG. 5: (Color online) The standard deviation A c of the 
charge on the circle k — 6 around the ring, as a measure of the 
amplitude of the charge deformation. Shown are the results 
for the ground state (GS), the 1-st excited states (ESI), and 
the 2-nd excited state (ES2), for TV = 2, 3, 4 electrons with- 
out without (tc = 0), and with interaction (tc = 2.2). The 
amplitude of the CDD's is strongly reduced by the Coulomb 
effects for TV = 3, 4; notice the different scales used in the 
paired panels (c),(d) and (e),(f). The magnetic energy ts 
varies between and 10 and the lines corresponding to the 
excited states are shifted to the right, for clarity. 



only generates the specific effective potential which de- 
termines the CDD. In the absence of SOI (a = (3 = 0), 
we checked that a similar screening effect occurs in the 
presence of a potential that induces a charge deformation 
comparable to that obtained with the SOI. It is, however, 
surprising that by adding only one extra electron such a 
drastic effect ensues. 

Next, we investigate the effect of the Coulomb inter- 
action on the spin distribution around the ring. The 
standard deviation of the spin density projected along 
the z direction, A z , is plotted in Fig. [6] where we show 
the results calculated as before for the circle correspond- 
ing to the sites with radial coordinate k = 6. The spin 
density deformation (SDD) is actually amplified by the 
Coulomb interaction for all TV = 2,3,4. As the CDD's, 
the SDD's reach their maximum at the magnetic fields 
where level repulsion occurs and remains prominent even 
when the gaps are very small. The Coulomb enhance- 
ment is a result of the mixing of states with the same 
parity, but with different spin orientation produced by 
the Coulomb potential. Consequently the SDD's have in 
general a richer structure than the CDD's. 

Finally, in Fig. [7] we display an example of CDD and 
SDD, obtained for TV = 4 interacting particles. The 
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charge and spin distributions corresponding to the sec- 
ond excited state and tg = 4.5 are illustrated in Figs. 
[5jf) and[6]Jf), respectively. The CDD is weak, but still it 
has four visible maxima. For two electrons the CDD has 
only two maxima which are along the directions x = y 
or x = —y, depending on the state and on the magnetic 
field, both with and without the Coulomb interaction. In 
particular, for a strictly one-dimensional ring model and 
N = 2, in the ground state, the maxima are always along 
the line x = —y^^ whereas for a two-dimensional model 
they can also be along x = y~ But in general, for N > 2 
electrons, screening effects may distribute the charge in 
more complicated configurations. Similar profiles with 
multiple local oscillation may be obtained for the spin 
density, eventually becoming spin-density waves around 
the ring. 
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FIG. 6: (Color online) The standard deviation A z of the spin 
projection along the z direction on the circle k — 6 around 
the ring, as a measure of the amplitude of the spin-density 
wave. Like in the previous figures GS, ESI, and ES2 in the 
legend indicate the ground state, the 1-st excited states, and 
the 2-nd excited states, respectively. The results are shown for 
N — 2,3,4 electrons without {to = 0), and with interaction 
(t c = 2.2). Unlike the CDD's, the SDD's are amplified by 
the Coulomb interactions for all N. The magnetic energy ts 
varies between and 10 and the lines corresponding to the 
excited states are shifted to the right, for clarity. 
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FIG. 7: (Color online) (a) The charge density for iV = 4 
electrons with interaction (tc = 2.2), in the second excited 
state, i. e. ES2 in FigOf), with magnetic energy energy ts — 
4.5. (b) The corresponding total spin distribution along the 
ring k — 6 where the standard deviation of the z component 
is calculated and shown in FigJBff). 



IV. CONCLUSIONS 



We calculated the many-body states of a system of 
N = 2,3, and 4 interacting electrons located in a ring of 
finite width with Rashba and Dresselhaus spin-orbit cou- 
pling, in the presence of a magnetic field perpendicular on 
the surface of the ring. The Coulomb effects are fully in- 
cluded in the calculation via the "exact diagonalization" 
method. Wc obtained inhomogeneous charge densities, 
or CDD's, around the ring due to the combined effect of 
the two types of SOI. When the Coulomb interaction is 
included the charge deformation is amplified for N = 2, 
as also shown by other authors For N > 2 we find 
that the CDD is dramatically flattened out in the pres- 
ence of the Coulomb interaction. We interpret the result 
as a screening effect. On the contrary, the spin inhomo- 
geneities, or SDD's, are amplified by Coulomb effects for 
all N > 1. 



Acknowledgments 

This work was supported by the Icelandic Research 
Fund. Valuable discussions with Sigurdur Erlingsson, 
Gunnar Thorgilsson, and Marian Ni^a are cordially ac- 
knowledged. 



1 S. Datta and B. Das, Appl. Phys. Lett. 56 (1990). 



2 Y. Bychkov and E. I. Rashba, JETP Lett. 39, 78 (1984). 



7 



3 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

4 R. Winkler, Spin orbit coupling effects in two-dimensional 
electron and hole systems ( Springer- Verlag Berlin, Heidel- 
berg, New York, 2003). 

5 T. Ihn, Semiconductor nanostructures. Quantum states 
and electronic transport (Oxford University Press, 2010). 

6 J. Splettstoesser, M. Governale, and U. Ziilicke, Phys. Rev. 
B 68, 165341 (2003). 

7 S. Souma and B. K. Nikolic, Phys. Rev. B 70, 195346 
(2004). 

8 J. S. Sheng and K. Chang, Phys. Rev. B 74, 235315 (2006). 

9 Y. Liu, F. Cheng, X. J. Li, F. M. Peeters, and K. Chang, 
Phys. Rev. B 82, 045312 (2010). 

M. P. Nowak and B. Szafran, Phys. Rev. B 80, 195319 
(2009). 

1 F. E. Meijer, A. F. Morpurgo, and T. M. Klapwijk, Phys. 
Rev. B 66, 033107 (2002). 



The first and second derivatives of any 
generic function f(x), are approximated as 
f'{x) « [f(x + h)-f(x-h)]/h and f"{x) » 
[f(x + h) + f{x — h) — 2f(x)] /h 2 , respectively, where 
h is considered sufficiently small. 

3 P. Hawrylak and D. Pfannkuche, Phys. Rev. Lett. 70, 485 
(1993). ' 

4 N. T. T. Nguyen and F. M. Peeters, Phys. Rev. B 83, 
075419 (2011). 

5 N. T. T. Nguyen and F. M. Peeters, Phys. Rev. B 78, 
045321 (2008). 

6 D. Pfannkuche, V. Gudmundsson, and P. Maksym, Phys. 
Rev. B 47, 2244 (1993). 

7 J. Schliemann, J. C. Egues, and D. Loss, Phys. Rev. Lett. 
90, 146801 (2003). 



